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چکیده 


چگونگی تغیرات غلظت رسوب معلق در طول مسیر رودخانه و بررسی عوامل تأثیرگذار بر آن همواره مورد توجه مهندسین علم هیدرولیک و محیط 
زیست می‌باشد. عدم وجود ایستگاه‌های اندازه‌گیری کافی و مشکلات نمونه‌برداری رسوب معلق» تهیه مدل‌هایی که به درستی رسوب معلق را در طول 
مسیر رودخانه روندیابی نمایند ضروری می‌نماید. در طبیعت رسوب بصورت نامتعادل انتقال می‌بابد در حالی که خیلی مدل‌های تجاری حالت ظرفیت حمل 
یا متعادل را در نظر می‌گیرند. از این رو در تحقیق حاضر مدلی عددی تهیه شده با حل عددی معادله تک‌بعدی انتقال و پخش غیرماندگار روندیابی 
رسوب معلق در یک بازه رودخانه‌ای را در شرایط عدم تعادل انجام می‌دهد. پس از صحت‌سنجی مدل, تأثیر ده روش عددی منفصل‌سازی» پنج معادله 
انتقال رسوب» هشت رابطه ضریب پخشیدگی و هشت رابطه سرعت سقوط ذره بر تغییرات بار رسوب معلق در طول بازه مورد مطالعه بررسی شد. نتایج 
تحقیق نشان داد استفاده از رابطه تحربی وایف مقدار رسوب معلق بیشتری را نسبت به سایر روابط دیگر برآورد می‌کند. در میان روابط سرعت سقوط ذره 
رابطه استوکس سرعت سقوط بیشتری را برآورد می‌کند که باعث می‌شود احتمال معلق شدن ذرات رسوب کمتر و در نتیجه غلظت آن نسبت به سایر 
روش‌ها کمتر باشد. همچنین در بین روش‌های منفصل‌سازی روش وان لیبر خطای کمتری را دارا است. از طرفی رابطه الدر کمترین و رابطه کاشفی 
پور- فالکونر پیشترین مقدار پخشیدگی را در هیدروگراف غلظت از خود بجا می‌گذارند. علاوه بر این نتایج تحقیق حاضر نشان داد غلظت رسوب برآورد 
شده در حالت عدم تعادل حدود ۱۱/۷ درصد بیشتر از ظرفیت حمل بار معلق محاسبه شده توسط روابط تجربی می‌باشد. 


واژه‌های کلیدی: روندیابی رسوب معلق. سرعت سقوط شبیه‌سازی عددی. ضریب پخشیدگی, معادله انتقال و پخش 


مقدمه 


امروزه یکی از مهمترین نگرانی‌ها در مسائل موجود در رودخانه‌ها 
عوامل تاثیرگذار بر میزان رسوب ورودی به مخازن سدها و دریاچه‌ها 
بارمعلق ورودی می‌باشد. برآورد دقیق میزان بار رسوب معلق از 
جنبه‌های گوناگون اجتماعی و اقتصادی حائز اهمیت است. 

با توجه به اهمیت موضوع تاکنون روش‌های مختلفی برای 
محاسبه بار رسوب معلق ارائه شده است. این روش‌ها می‌توانند مبنای 
توفالتی (۱۰) و وان‌راین (۱۱)). این معادلات ارائه شده عموما ظرفیست 


۱ و ۲- به‌ترتیب دانشیار و دانشجوی کارشناسی ارشد سازه‌های آبی» گروه مهندسی 
آب. دانشگاه رزی» کرمانشاه 

جرمم, تمصع ۵ مفلدها0طع15:اتقصصط) 
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(#- نویسنده مسئول 


حمل رسوب به صورت بار معلق را برآورد می‌نمایند. ژانگ و همکاران 
(۱۳) اعلام داشتند بسیاری از روابط انتقال رسوب برای شرایط تعادلی 
توسعه داده شده‌اند که در این حالت فرض بر این است که ظرفیت 
حمل رسوب با مقدار انتقال رسوب برابر است. با این وجوده وقتی که 
بستر کانال در حال فرسایش و رسوب‌گذاری باشد. میزان انتقال رسوب 
ممکن است کمتر پا بیشتر از ظرفیت حمل رسوب گردد. 

به منظور برآورد رسوب معلق در شرایط عدم تعادل که در آن 
تباال ذرات رسوب بین بستر و ستون آب با تاخیر مکانی و زمانی اتفاق 
می‌افتد و به صورت آنی نیست. لازم است در معادله اتتقال بای ۲ 
کف زانط وی را میاه غلظه فیاکین شطع 
همراه با ترم منبع به صورت عددی حل شود. 
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پسته با شبیه‌سازی نوع آلودگی. ترم منبع 80) از روابط مختلفی 
محاسبه می‌شود. که برای رسوب معلق می‌توان از رابطه زیر استفاده 


کرد. 
۳ 3 ِ 7 8 
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که در آن 0 غلظت رسوب معلق» ,۵ سرعت سقوط ذرات 
رسوب. ,7 عرض سطح آب؛ ,4 ظرفیت حمل رسوب معلق در 
واحد پهنای کانال» ,4 دبی جریان در واحد عرض کانال 1 ضریب 
وان کارمن که معمولا برابر ۰/۴۱ در نظر گرفته می‌شود و *[) 
سرعت برشی می‌باشد. حل معادله ۱ تحقیقات زیادی را به خود 
معطوف داشته است از جمله آپادو (۱) که با استفاده از سه روش لکس 
وندروف به صورت صریح و کرنک نیکلسون به صورت ضمنی و روش 
تفاضل محدود غیر استاندارد به حل معادله یک بعدی جابجایی و 
پخش با ضرایب ثابت پرداخت. نتایج ایشان نشان داد روش‌های لکس 
وندروف و روش تفاضل محدود روش‌های مناسبی برای حل معادله 
یک‌بعدی جابجایی و پخش است. 

کاشفی‌پور و توکلی‌زاده (۴) به معرفی مدل یک بعدی 
هیدرودینامیکی و کیفی فستر" و همچنین استفاده از روش کرنک 
نیکلسون به شکل تفاضل مرکزی پرداختند. پس از اجرای مدل برای 
بکتصال عامن مت یرای انفان ماکان 
۵ و مقایسه با مقادیر تحلیلی بدست آمده مدل فستر را مدلی دقیق 
اعلام کردند. باغبان‌زاده و کاشفی‌پور (۲) به کمک مدل عددی و 
استفاده از معادلات سنت-ونانت و هم چنین معادله انتقال-پخش 
برای مدل‌سازی جربان و انتقال رسوب معلق به مطالعه انتقال رسوب 
معلق غیر چسبنده در رودخانه کرخه پرداختند. پس از اجرای برنامه با ۵ 
معادله رسوب معلق و ۶ معادله برای تعیین ضریب پخش و مقایسه با 
داده‌های اندازه‌گیری شده ترکیب معادله رسوب معلق وان‌راین و 
معادله ضریب پخش فیشر نتایج مناسبی را ارائه کرد. صدائی و 
همکاران (۷) به بررسی میزان دقت شش فرمول بگنولد لین و 
کالینسکی انیشتین» چانگ-سایموز-ربچاردسون, بروکس, توفالتی در 
برآورد بار معلق رسوب نسبت به داده‌های اندازه‌گیری شده در رودخانه 
سولگان پرداختند. با انجام این تحقیق مشسخص شد که از بین 
فرمول‌های بررسی شده فرمول بروکس برای سولگان دارای برآورد با 
دقت بالایی بود. صدائی و سلیمانی (۸) به بررسی دو فرمول لین و 
کالینسکی و بگنولد برای محاسبه بار معلق در رودخانه ارمند با مقادیر 
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اندازه‌گیری شده و اجرای روش تحلیل مسیر در رودخانه ارمند 
پرداختند. در نهایت با بررسی‌های صورت گرفته مشسخص شد که 
فرمول بگنولد دارای برآورد نزدیکتری نسبت به فرمول لین و 
کالینسکی است. صداتی و سلیمانی () به ارزیابی فرمول‌های برآورد 
بار رسوب معلق (توفالتی» چانگ-سایمونز -ریچاردسون و بگنولد) بر 
مبنای داده‌های سه رودخانه بازفت» سولگان و ارمند با استفاده از 
روش تحلیل سلسله مراتسی" پرداختند. فرمول چانگ-سایمونز- 
ریچاردسون در رودخانه‌های دکر شده به عنوان فرمول مناسب جهت 
برآورد بار رسوب معلق پیشنهاد گردید. همچنین حبیبی (۳) بر اساس 
مفاهیم اصلی هیدرولیک انتقال انرژی و حرکت ذرات در جریان آب و 
رسوب و با فرض اینکه غلظت رسوبات در هر عمق مشخص از کف 
کان است نا ری شیاه شنم در ماد سم 
می‌باشد یک رابطه نیمه تجربی ارائه داد. 

لازم به ذکر است در سال‌های اخیر تحقیقات زیادی مبنی بر 
هوش محاسباتی غلظت بار معلق در مجاری روباز از جمله پژوهش 
رجایی و میر باقری (۶) که با استفاده از شبکه عصبی مصنوعی به 
پیش‌بینی غلظت بار معلق پرداختند انجام گرفته است. پس از مقایسه 
با روش منحنی سنجه رسوب. نتایج نشان‌دهنده برتری مدل‌های 
شبکه عصبی مصنوعی بود. 

مطالعات قبلی عمدتا بر روی تعیین نسبی مقدار بار مق و 
مقایسه آن با روابط ارائه شده مبنی بر ظرفیت حمل در محل یک یا 
چند ایستگاه هیدرومتری در یک رودخانه اعلام شده است. بررسی 
تغییرات غلظت بار معلق و عوامل موثر بر آن در طول یک بازه‌ای از 
رودخانه به ویژه در شرایط عدم تعادل که به واقعیت نزدیک‌تر است در 
این تحقیق مورد توجه قرار گرفته است. از این رو در تحقیق حاضر 
بازه ورودی رودخانه قره‌سو به شهر کرمانشاه که از نقطه‌نظر محیط 
زیستی بسیار حائز اهمیت است مورد توجه قرار گرفت. با این وجود 
به علت نبود ایستگاه اندازه‌گیری رسوب و جریان در محدوده مذکور و 
همچنین مشکلات و هزینه بسیار زیاد اندازه‌گیری بار رسوب معلق از 
طرف دیگر» لزوم تحقیقی به منظور بررسی تاثیر پارامترهای مختلف 
بر توزیع غلظت بار معلق در طول بازه مورد نظر و تعیین محدوده‌های 
بحرانی از نقطه نظر آلودگی رسوب معلق را دو چندان می‌نماید. 


مواد و روش‌ها 
محدوده مورد مطالعه 

برای مدل‌سازی رسوب معلق در اين مقاله. رودخانه قره‌سو به 
طول ۱۱۴۸۰ متر و تعداد ۰ مقطع عرضی برداشت شده مورد بررسی 
قرار گرفته است. همانگونه که در شکل ۱ نشان داده شده است جهت 
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بررسی عددی عوامل موّثر بر توزیع غلظت رسوب معلق نامتعادل در رودخانه‌های طبیعی ۲۳۳ 


کلی جریان در سیستم رودخانه‌ای مذکور از شمال غرب به سمت 

جنوب شرق می‌باشد. برخی مقاطع عرضی برداشت شده در شکل ۲ 
نشان داده شده است. در رودخانه قره‌سو چهار دانه‌بندی مواد بستر 
برای بازه مورد مطالعه در نظر گرفته شده است. در شکل ۲ منحنی 
دانه‌بندی مواد بستر ارائه شده است. در این شکل دانه‌بندی تیپ ۱ 
برای مقطع عرضی ۱ تا ۱۰ (کیلومتر صفر تا ۲۷۵۴ دانه‌بندی تیپ ۲ 
برای مقاطع ۱۱ تا ۲۶ (کیلومتر ۲۷۵۴ تا ۶۵۰۵) دانه‌بندی تیپ ۳ 
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دانه‌بندی تیپ چهار برای مقاطع ۴۲ تا ۶۰ (کیلومتر ۱۲۱۲۶ تا 
۵ استفاده شده است. از دانه‌بندی و مقاطع عرضی به عنوان 
شرایط اولیه ورودی مدل جریان و رسوب استفاده شد. 

همچنین تغییرات ضریب زبری در طول بازه, مطابق با جدول ۱ 
برای ساحل چپ. کانال اصلی و ساحل راست رودخانه برای مدل در 
نظر گرفته شده است. 
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شکل ۱- پلان رودخانه قره‌سو 
۱۱۷۵ ممقحتقطای ۵۶ »تافصمد 5 1 مساع۲1 
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شکل ۲- پروفیل عرضی برخی از مقاطع رودخانه قره‌سو 
۵۲ 22600۵ ۵۴ دممتا0ع5 عومن 506 2 ۲12۱0۲6 


جدول ۱- مقادیر مربوط به ضریب زبری مانینگ رودخانه قره‌سو در بازه مورد مطالعه 
2 5۲0۱0۱۵4 6ظ۳) مد جهن ومعه دای ۵ هی وعمصطامت0 مسنصطه]۱۷ -1 ماماو] 


فاصله از پایین‌دست (متر) ضریب ساحل چپ ضریب کانال میانی ضریب ساحل راست 
[۳] همهم صم معصه) مزر اقصه ]1,۵[ 2 امه ۱۷۲۵10 2 اقصی )1:۵۳ ۵ 

0045 036 0045 3396.639 0 00 
0045 0030 0045 3396.639 ۱0 7 
0.050 0028 0050 6374857 0 03 
0.040 ۱0۳0۶۸ 0040 10520.45 0 8 
00048 0030 00048 11714:28 0 1 
0045 025 0045 13600.91 0 3 
0045 0030 0045 15026.23 6 0 
0045 025 0045 15909.50 16 3 
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شکل ۳- دانه‌بندی مواد بستر در طول رودخانه قره‌سو 
۵۲ 621219500 هصماه ع۲ین متا طاتتاعنل 5126 نی -3 م۲۵ 


معادلات جریان 

اولین مرحله در محاسبات بار معق تعیین خصوصیات 
هیدرولیکی جریان در مقاطع مختلف می‌باشد. از آتجائی که در این 
تحقیق شرایط جریان ماندگار و شبه ماندگار تک‌بعدی در نظر گرفته 
شده است برای محاسبه نیمرخ سطح آب در اين رایط از روش گام 
به گام استاندارد و حل رابطه انرژی برای مقاطع نامنظم به عنوان 
معادله حاکم استفاده شد: 

2 2 

(۳/ وگو وه ۷ و < 2 ج و ره +۱۱ 

۸ سطح مقطع. 0 دبی جریان, 2 رقوم کف از عمق جریان 
و 67 میانگین شیب خط انرژی است که برابر: 


۵( ت رک 7 
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- و 
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که درآن 7 ضریب زبری مانینگ و1 شعاع هیدرولیکی 
مقطع شماره ۲ در پایین دست با فاصله بر از یکدیگر دارد. در 
تحلیل مسائل در کانال‌های با سطلح مقطع منظم 1 :6 فرض 


محاسبه گردید. 
2 3 
0 ‌ 2 ۳ - 
(ب4ع) 
1 - ۸ 
)۵ تب 


در این روابط ,/ ضریب انتقال در زیر مقطع / ام از مقطع مورد 
نظر می‌باشد. 

پس از تغییر رابطه ۴ به شکل رابطه ۱۰ بامعلوم بودن 
مشخصات جریان در مقطع کنترل (مثلا مقطع شماره ۱) و همچنین 
وارد کردن مقاطع عرضی به مدل با استفاده از روش نصف کردن 
فاصله و حل رابطه غیرخطی زیر بر حسب ۷ مقادیر عمق جریان 
و در نهایت سرعت و نیمرخ طولی سطح آب در طول بازه محاسبه 

0۶ 
۳ 2 * رو ,۵ ,۷ < (ومر)۶ 
2 


0 - ۸ 5۴ + و2 + چپروچ 0 + وم 
2 
معادلات رسوب 


با استفاده از روش احجام محدود منفصل‌سازی معادله ۱ به 
روش کاملا ضمنی انجام شد و رابطه جبری زیر بدست آمد: 


0 
+ + 0 + 
(«م باه ۵,0 +(01 0ام۸ 
۵ 0۳ 0۳7 0 0۳ 0 
520 ِ 0 775" ۳ كِ 5" ما۵ 


به منظور تقریب زدن مقادیر غلظت بر روی سطوح حجم 
کنترل ( ۳,0۳ ) در رابطه ۱۱ طرح‌های مختلفی وجود دارد. 
در این تحقیق از روش‌های توان -پیروء کوییک» اختلاف مرکزی و 
طرح کاهش اختلاف مجموع ‏ استفاده شده است. 
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جدول ۲- روش‌های منفصل‌سازی, معادلات سرعت سقوط. ضریب پخشیدگی و رسوب معلق 
جمناجنوی )صمصنلعی 060صعمرویاه_ 0ص )صعمن‌گلهم صمتمیگنظ۱ , مصمتاهاعد وان۷۵۱۵ 1211 وول ماع م6 هعتاع جوز( -2 12016 


۱ فرمول‌های سرعت 
فرمول‌های بار معلق فف ان 
سقو 

ص۲0 و 02 ۵ 
(1941) ع1وصلق! 6ص[ (1851) 506016 
(1963) 101016101 (1993) 0۷00۷ 
(1984) حززر و2 ۱۷ (1977) 7001 
(21.)1966 61 ۱۷۷۱۲6 (1984) عصهل 


(1984) ۲۱20101 (1989) طرز۲ ۲۷2 


(1995) صعتانا[ 
(1999) ب«هطاواناه5 
(1989) ع72 


علاوه بر روش‌های منفصل‌سازی معادله انتقال و پخشیدگی» 
روابطی که برای تخمین سرعت سقوط ذرات رسوب» ضریب 
پخشیدگی و روابط تخمین رسوب معلق مورد استفاده شده در این 
تحقیق به اختصار در جدول ۲ آورده شده است. 


شرایط مرزی 

مقاطع هندسی بازه مورد نظر و منحنی دانه‌بندی رسوبات به 
عنوان شرایط اولیه به مدل وارد می‌شود. شرط مرزی بالادست برای 
مدل جریان, مقدار دبی مشخص و شرط مرزی پایین‌دست رابطه 
دبی- اشل است که برای مقطع آخر رودخانه قره‌سو با توجه به شیب 
عمومی رودخانه از رابطه مانینگ ساخته می‌شود. 

در مدل انتقال-پخش شرط مرزی بالادست مقدار غلظت معلوم 
است این غلظت می‌تواند بصورت پیوسته آنی با مدت محدود از 
بالادست وارد رودخانه مذکور شود. شرط مرزی گرادیان صفر به 
عنوان شرط مرزی پایین‌دست برای مقطع آخر رودخانه قره‌سو تعریف 
شده است. مدل مورد توجه این تحقیق در محیط ویزوال بیسیک 
تشه که مازلا رتوب»صر مخ ظوان بوعت فرط 
و روش‌های منفصل‌سازی مختلف در زیر برنامه‌های این مدل تعریف 
شده که در هنگام اجرای برنامه بر حسب مورد فراخوانی می‌شوند. 


صحت سنجی مدل هیدرولیکی 

به منظور بررسی دقت مدل هیدرولیکی از نقطه‌نظر برنامه‌نویسی, 
از محدوده مورد مطالعه این تحقیق یعنی رودخانه قره‌سو با طول 
۲۰ متر و تعداد۶۰ مقطع عرضی استفاده شد. بدین منظور 
مشخصات رودخانه به همراه دبی متوسط ورودی ۲ ساله 2/5 
۴ به عنوان شرط مرزی بالادست و رقوم سطح آب که از رابطه 


ضریب پخشید گی روش‌های منفصل سازی 


اصمن‌نگلهی موز 1۲16۲6۱6۵6( 


(1974) ۵۵۲6۵۲ 220 ۷۲۵0۵1۷۵۷[ ۳0۷۵1۳ 
(1966) ۲1۱50۳6۲ 104 
(1959) ۲۱067 0 
(۳۵۱6۵۴۵۲)2002 20 کالوک 6۲ صه ۲۷ 
(1(112)2009 240 *2[66؟] 2 2 ۱۷ 
(1998) .21 )6 1[ 6 ۷۲1۴[ 
(1998) 06060782 200 50 ۱ 1۳/9 
(1967) 1[ 0 
0۵00 
اعتصرنا 


میک :هه می شوه یه عوان فرط مر زین این دلست بته:ملال 
حاضر و نرم‌افزار هکرس وارد شد. 

به منظور بررسی عملکرد مدل» هر دو مدل به [زای ضریب زبری 
تک کی رل وس نوی تسط آن یبای 
محاسبه شد. همانطور که در شکل ۴ مشسخص است نیمرخ طولی 
سطح آب محاسبه شده در دو مدل بر هم منطبق شده است به 
ری ک یگیم ای مقلی فرع ستلم آب شدانید تال 
حاضر از مدل هکرس حدود ۵/۲ سانتی‌متر در مقطم شماره ۱۱ و 
ميانگین خطای مطق ۱/۵۲ سانتیمتر در طول بزهمیبشد. 


صحت‌سنجی مدل رسوب 

به منظور بررسی دقت مدل شبیه‌ساز رسوب از نقطه‌نظر 
برنامه نویسیء با حفظ همان شرایط جریان اشاره شده در قبل. 
هیدرو گراف غلظت نشان داده شده درشکل ۵ به عنوان ورودی به 
مدل معرفی شد. قله اول هیدروگراف مذکور شبیه تزریق پیوسته و قله 
آخر شبیه تزریق لحظه‌ای غلظت در رودخانه است علاوه بر این قله 
وسط شبیه تزریق به مدت محدود است. حجم زیر این هیدروگراف 
نمودار غلظت-زمان ۱۸۲/۴۷۵ می‌باشد. با حذف ترم منبع در رابطه ۱ 
مدل به ازای روش‌های مختلف منفصل‌سازی اجرا شد. حجم زیر 
هیدروگراف‌ها دریک مقطع میانی بازه (مقطع ۱۵) و انتهای بازه 
(مقطع )۵٩‏ با ورودی مقایسه شد. به منظور مقایسه دقیق‌تر مدل 
مقادیر حجم زیر هیدروگراف‌ه ا به |زای روش‌های مختلف 
منفصل‌سازی در جدول ۲ ارائه شده است. ماکزیمم انحراف حجم زیر 
هیدروگراف مقاطع میانی و خروجی بازه از هیدروگراف ورودی ۰/۰۲۹ 
درصد است. اين مقدار کم خطا نان از دقت بالای مدل و اطمینان 
به نتایج آن برای ادامه تحقیق دارد. 
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ارتفاع (متر) 


ز[ض] صمنادبه؟ 


فاصله (کیلومتر) 
[۴2] ععصهاوزر[ 
تراز سطح آب محاسبه شده با مدل حاضر تراز سطح آب محاسبه شده با مدل هکرس 6 تراز بستر کانال اج ام 
,1 هناد تماه: ۲۵06 اهعوع‌رعط 0 0 آ2 1 ۱۶ .0 0عظ 


شکل ع- مقایسه پروفیل سطح آب محاسبه شده با مدل هکرس و مدل حاضر 
۵061 0۲66۵0۵ ۵00 ۲۱۳۰۸۸۹ ۱ ممتاهماه هه ماه 0 مفزهمج0) -4 عاع۲1 


10 
5 
ِِ 
2 
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0 
60 50 40 30 20 10 0 
زمان (ساعت) 
[ط] 11۳6 


شکل ۵- هیدروگراف غلظت ورودی به بازه مورد مطالعه 
۵ 5000۱60 ۲6 م6 حامره عم و متاحماصهعصم )سامرصا -5 مساهز۲ 


جدول ۳- حجم زیر نمودار غلظت_-زمان برای روش‌های مختلف در مقاطع ۱۵ و ۵٩‏ 
9 200 15 کصمناع6ز ص د0مط)مهه اصمعمن 0 جارهعم0 مطصمتاهصاصم صم معط مامت ممصتاه۲ 1 -3 12016 


مقطع ۵٩‏ مقطع ۱۵ 
9 «ما0ع 15 6007 
حجم زیر نمودار میزان خطا حجم زیر نمودار روش منفصل‌سازی 
حاصرحه00 ۲ ماه تامیمصت مصصسان ۱۲۷ ۵01۳0 ۹ ۰ وه با 

199240 00018 1999240 ۱۹ 
18240 0018 18240 4 
18240 0018 18240 امن 
19251 0029 19251 ۲ 22 ۱۷ 
18240 +0024 182.50 2 2 ۲۷ 
18240 0029 19251 ۷۲10-06[ 
182.50 0۳94 182.50 0 ما5 
18240 0013 199248 7 
18240 0018 18240 ۹ 
18240 +0024 182.50 اعتصرل] 
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احناجعن) 0۷۷۵۲ 
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1 1 1 0 ۹ ۰ ۰ 0 
2 
60 20 0 60 40 20 0 0 ب وا 20 9 
ان ات زمان(ساعت) 
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شکل ۷- هیدروگراف غلظت محاسبه شده به روش‌های مختلف منفصل‌سازی 
عمصامصه حمتاهمتاع معا )مهن وا حامرصهم بط ممتاح‌اصصمعصی 0ع)علنهامن -6 ۲1۵0۲۵ 


تعیین مناسبترین روش منفصل‌سازی داده شده در شکل ۵ (که به عنوان ورودی در این حالت نیز در نظر 
برای تعبین بهترین روش منفصل‌سازی, با حذف ترم منبع و ترم گرفته می‌شود) تنها در جهت محور زمانی جایجا نماید و هیچگونه 


پخشیدگی و پایین افتادگی در قله‌های هیدور گراف رخ ندهد. به منظور 
بررسی دقیقتر برای روش‌های مختلف منفصل‌سازی اشاره شده در 
جدول ۲ مدل اجرا گردید. همانگونه که در شکل ۶ مشاهده می‌شود 
همه روش‌های منفصل‌سازی به جز طرح اختلاف مرکزی» برای 


پخشیدگی» رابطه ۱ به رابطه انتقال خالص" تبدیل می‌شود. بهترین 
روش منفصل‌سازی روشی است که در آن هیدره گراف غلظت نشان 


1- 0۷1۳۶ 8 
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شبیه‌سازی قله اول و دوم هیدروگراف ورودی از دقت مناسبی 
برخوردار هستند ولی در شبیه‌سازی قله سوم هیدروگراف که شکل 
نوک تیز دارد و رعایت الگوی هیدروگراف ورودی روش وان لییر از 
دقت بیشتری برخوردار است. 


پس از صحت‌سنجی مدل‌های جریان و رسوب و اطمینان به 
نتایج آنها تأثیر ضریب پخشیدگی بر توزیع غلظت بررسی شد. بدین 
منظور مجموعه‌ای از روابط ارائه شده در منابع مختلف برای تعیین 
ضریب پخشیدگی جمع‌آوری شد و با معرفی هر یک از آنها به مدل, 
هیدروگراف غلظت در خروجی بازه محاسبه گردید. همانگونه که در 
شکل ۷ نشان داده شده است استفاده از رابطه الدر برای محاسبه 
ضریب پخشیدگی و جایگزینی در معادله انتقال و پخش [رابطه ۱) 
باعث شده است که قله اول» دوم و سوم هیدروگراف کمترین اختلاف 
غلظت و تغییر شکل را با هیدروگراف وروی به بازه داشته باشند. به 
عبارتی رفتار معادله انتقال و پخش بیشتر شبیه معادله انتقال خاص 
شده است و دلیل این امر مقادیر ناچیز ضریب پخش محاسبه شده 
توسط رابطه الدر می‌باشد. از طرفی رابطه کاشفی‌پور و فالکونر 
بیشترین مقدار پخشیدگی را در هیدروگراف غلظت از خود بجا گذاشته 
است. استفاده از این رابطه باعث شده است که هر سه قله هیدروگراف 
غلظت فروافتادگی داشته باشد و به دلیل حفط پیوستگی (حجم زیر 
هیدرو گراف‌های غلظت در خروجی با هیدروگراف ورودی برابر است) 
هی هیفروگ اف اهنا بان بانط کی را ورد 


است. 


مقایسه روابط سرعت سقوط 


سرعت سقوط ذره یکی از پارامترهای تاثیرگذار در انتقال رسوب 
می‌باشد. از آنجائیکه روابط متعددی برای محاسبه سرعت سقوط ارائه 


ده وه نامه( 
۱ 
هل جح 
2160 ۳ مج عبرم حرتله دوه 1 
۱۹ 
6 1 
60۶ 220 560 
۳ 


40 50 0 


30 
زمان(ساعت) 
[ط] م1 
شکل ۷- تأثیر ضرایب مختلف پخشیدگی بر توزیع غلظت رسوب 
مماطا تنل )اعد جم کاصهت‌الهی صمتعدگلنل ]0۱1 0۲ ۲۱2۱۷۲۵-۲۳۲6 


جمع‌آوری و تأثیر آنها بر غلظت رسوب رودخانه قره‌سو مورد بررسی 
قرار گرفت. بدین منظور از روش منفصل‌سازی وان لییر به عنوان 
روش برتر رابطه مک کوثیری و کیفر برای ضریب پخش طولی و 
همچنین استفاده از رابطه بار معلق وان‌راین برای محاسبه ظرفیت 
حمل رسوب معلق استفاده شد. در شکل ۸ توزیع غلظت رسوب بار 
معلق در طول بازه به ازای روابط مختلف محاسبه سرعت سقوط نشان 
داده شده است. لازم به ذکر است با توجه به غیریکنواختی دانه‌بندی 
رسوب در طول بازه غلظت رسوب بار معلق برای هر ذره رسوب 
جداگانه محاسبه شده و سپس با توجه به درصد وزنی آن ذره در 
منحنی دانه‌بندی غلظت رسوب کل ذرات( بصورت حجمی و بدون 
بعد) در هر مقطع محاسبه و در شکل ۸ نشان داده شده است. 


همان‌طور که در شکل مشخص است کمترین غلظت رسوب بار 
معلق هنگامی بدست می‌آید که از رابطه استوکس برای محاسبه 
سرعت سقوط استفاده شده است. استفاده از رابطه استوکس سرعت 
سقوط بیشتری در مقایسه با سایر روش‌ها ارائه می‌دهد و احتمال 
معلق شدن ذرات را کمتر نشان می‌دهد. روابط وان‌راین و زانک در 
مقایسه با سایر روابط غلظت رسوب بیشتری را نشان می‌دهد. 


تأثیر معادله رسوب معلق 

به منظور بررسی تاثیر معااله رسوب معلق در برآورد غلظت 
رسوب در طول بازه مورد بررسی؛ برخی روابط تجربی مشهور به 
عنوان ,4 در رابطه ۲ در نظر گرفته شد. بدین منظور از روش 
منفصل‌سازی وان لییر رابطه مک کوثیری و کیفر برای ضریب پخش 
طولی و سرعت سقوط روبی استفاده شد. تغییرات طولی غلظت 
متوسط رسوب در هر مقطع در شکل ٩‏ نشان داده شده است. 
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شکل ۸- تأثیر روابط مختلف سرعت سقوط بر توزیع غلظت 
ما متاعنل حمصصت0عد حدم فامم] چات۲۵۱۵ )0۱116۳6 0۴ ۳/۲۵۵۲ -8 ۲12۱۷۲۵ 
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شکل -٩‏ تأثیر روابط مختلف تجربی برآورد ظرفیت حمل رسوب بر توزیع غلظت 


ص تنل اصمصت0عی مه ی )نمی 060صمصرمیاه که دی )60ج و6 ممتاهاه همرت )هلت ۵۲ ۳/6۵۲ -9 مسوز۲ 


همانگونه که ملاحظه می‌شود روش‌های توفالتی» لین -کلینسکی» 
وان‌راین و حبیبی روند تقریبا یکسانی را نشان می‌دهند هرچند که 
روش وان‌راین مقدار بیشتری رسوب معلق را در عمده مسیر بجز 
حوالی کیلومتر ۳ در مقایسه با این سه روش برآورد می‌کند. علاوه بر 
این روش وایف بیشترین و تفالتی کمترین مقدار بار معلق رسوب را 
در بازه مورد مطالعه برآورد می‌کند. از جمله دلایل اختلاف فاحش 
روش وایف و همکاران با سایر روش‌ها در این است که در این روش 
ظرفیت حمل رسوب حساسیت بیشتری به سرعت جریان عمق و 
شیب رودخانه دارد به گونه‌ای که با توان دوم حاصلضرب این 
پارامترها در ارتباط است. 


مقایسه حالت تعادل و عدم تعادل 

غلظت رسوب معلق در حالت تعادل با استفاده از روابط تجربی 
ارائه شده توسط محققین مختلف برآورد می‌شود. غلظت رسوب در 
حالت عدم تعادل با استفاده از حل عددی معادله ۱۱ بدست می‌آید که 
می‌تواند بیشتر یا کمتر از حالت تعادل باشد. در شکل ۱۰ مقدار بار 
معلق را در حالت تعادل و عدم تعادل به ازای معادله‌های مختلف بار 


معلق نشان داده شده است. غلظت رسوب برآورد شده در حالت 
نامتعادل در مقایسه با ظرفیت حمل رسوب معلق محاسبه شده توسط 
روش‌های وایف و همکاران لین و کالینسکی, وان راین» حبیبی و 
تفالتی به ترتیب ۰۱۲/۷ ۰۱۱/۵ ۰۱۳/۹ ۱۲/۷ و ٩/۶‏ درصد بیشتر است. 

غلظت رسوب برآورد شده در حالت عدم‌تعادل همانگونه که در 
شکل ۱۱ نیز نشان داده شده است حدود ۱۱/۷ درصد بیشتر از ظرفیت 
حمل بار معلق می‌باشد. از آنجایی که در طبیعت؛ رسوب معلق بر 
اساس شرایط عدم‌تعادل منتقل می‌شوده عدم در نظر گرفتن این 
اختلاف در مسائل مرتبط با مدیریت رسوب از قبیل محاسبه عمر مفید 
سدها باعث می‌شود عمر مفید غیر واقع‌بینانه بیشتر برآورد شود. 


نتیجه گیری 

در این تحقیق مدلی عددی برای شبیه‌سازی انتقال رسوب معلق 
با دانه‌بندی غیریکنواخت و در شرایط عدم تعادل توسعه داده شد. 
صحت‌سنجی مدل در بخش هیدرولیک جریان و رسوب. نتایج دفیق 
و قابل قبولی را ارائه داد. در این مدل با حل عددی معادله انتقال و 
پخش غیراندگار همراه با ترم منبع» توزیع غلظت در طول بازه 


۰ نشریه آب و خاک. جلد ۳۴. شماره ۲. خرداد - تیر ۱۳۹۹ 


رودخانه قره‌سو محاسبه شد. در ترم منبع روابط تجربی متعددی برای 
برآورد ظرفیت حمل بار رسوب معلق در نظر گرفته شده است. ننایج 
یی شاخ ناه عاطت :نار رسوت شام :هر طول بانط هام فاد 
از روش‌های توفالتی» لین -کلینسکی» وان‌راین و حبیبی در ترم منبع 
تقریبا یکسان است. 

با این وجود روش وایف بیشترین و تفالتی مقدار کمترین غلظت 
رسوب معلق را برآورد می‌کند. تحقیقات قبلی نشان داده است که 
روابط انتقال رسوب مختلف با توجه شرایط هیدرولیکی» هندسی و 
برای هر رودخانه رابطه و يا روابط ویژه‌ای پیشنهاد شده است. برای 
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مثال باغبان زاده و کاشفی پور(۲) در شبیه سازی انتقال رسوب معلق 
غیرچسبنده در رودخانه کرخه مناسبترین رابطه شبیه‌سازی بار معلق را 
رابطه وان‌راین پیشنهاد کردند. صدائی و همکاران (۷) رابطه بروگز را 
برای شبیه‌سازی رسوب معلق رودخانه سولگان پیشنهاد کردند. صدایی 
و سلیمانی (۸) رابطه بگنود را در محاسبه رسوب معلق رودخانه ارمند 
انیب اور امه کیت و کشک کین کات تافو 
سلیمانی (*) با استفاده از روش تحلیل سلسله مراتبی فرمول برآورد بار 
رسوب معلق چانگ-سایمونز-ریچاردسون را برای سه رودخانه بازفت؛ 
سولگان و ارمند پیشنهاد کردند. 
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شکل ۱۰- مقایسه حالت تعادل و عدم‌تعادل با روابط مختلف رسوب معلق 
عصمناهنوع من )م0 وا مممتان0هصی صمترماناننمع-صمه مه انم ۵۶ صموتنهم‌صم -10 ع۲ع۲1 
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شکل ۱۱- مقایسه بار رسوب محاسبه شده در حالت تعادل و عدم تعادل 
کصمانصی صسترطناننوع-صمه مه ینانوی 0ص 0عم۱ 0ع0صهمرکیای اهامای ۵ جمونهم‌صصمن 11 عساع1ز۲ 


بررسی عددی عوامل موّثر بر توزیع غلظت رسوب معلق نامتعادل در رودخانه‌های طبیعی ۲۵۱ 
بدلیل اینکه سرعت سقوط محاسبه شده توسط رابطه استوکس از غلظت ایجاد کردند. باغبان زاده و کاشفی پور (۲) رابطه ضریب 
سایر روابط بیشتر برآورد گردید» غلظت رسوب معلق محاسبه شده در پخشیدگی فیشر را در شبیه‌سازی رودخانه کرخه پیشنهاد کردند. 
بازه هتگام در نظر گرفتن این رابطه از سایرروش‌ها کیتر بدست آمد. ۰ . بررسی میزان رسویات در حالت تعادل و عدم‌تمدل تشن داد که مقدار 
همچنین در بین روابط ضریب پخشیدگی طولی رابطه الار کمترین و رسوب برآورد شده در رودخانه قره سود در حالت عدم تعادل ۱۱/۷ 
رابطه کاشفی‌پور - فالکونر بیشترین مقدار پخشیدگی را در هیدروگراف درصد بیشتر از مقدار آن در حالت تعادل است. 
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۶ مدمه عمتهلنملهه 10۲ فممتاقام: 37 ۵۶ اععه راملمه مطا عصتهل۷20 عنم .ع0مطامصه ماه ها 00ومطهم 
عط فافع ممتلقناوه ۷۷1۲۵ مط عصوتا فطع 0ع«مطو فملناوع م1 ,لماهعلا9ع1۵۷ ۳/28 1۵20 )صمصصتلعو 0ع0هممفناه 
٩2) 0‏ 0وله ممتمناوه 101121610 م1 .فومتتهتاوع قطان صقطا تمطونط اصمحصتلعو 0ع4صوموفیه ۵۶ امه 
متطعصم‌تاقاه: وماماد رفحم؟ ات‌ماع۷ ام ملمتانهم )صعته عجمصه مصمتامتاوه تقطاه صقطا ها 1۵20 اصمحصت86 
عمط 0606۵865 صمئلهء‌ناموه صمتلدتاوع وعاماو فطا رععصق! .ععمتنمناوه قطانم صقط مومع لماع ۶ مطا 0متفصتاوم 
0 ۲۵۱2009919 معَلصعم2 مه منک ۷2 عصارمامره ره ۳۱۵۱۷۵۷ .عمامتاتهم اممصصتل0عو مطا هصت‌مممویاه ۵۶ وانلزماتوووم 
تنل 10۲ ومتطمممتاهام: اوه عممصصض ,جموع۲ تم فطع عمماه ممتاطتتاوتل 1۵20 عمممصتلعی 60صمم‌وناه تماقعتم 2 10 
۶ اصامصصح )فمطونط معط هه و10۳۵ مطا رصان فممتفنوع ممممله۳ - سبامم‌تلتطوفک مطا مه علاط رعامعته‌تااههه 
تممصنلعی ۵60همممناه ‏ 12060تملهه ‏ فص رمعمممطتا ۲‏ نموم هه ممتاحتاموعجمه _ فط صا صمتفتل 
عصمه عمتجم حصیت‌طاتاننوه تعلصه فطع صقطع ععطعنط 2۵ 11.7 ۷26 فممتات0جمع ترامع مه مهم ممتاه‌تامهم‌صمم 
6۰ 11۷۵۲ 36 

عص ]زمحفصها )حمحصنلعو وع1۵ لهاما مج ۵0 مه متقلیاصته ده واممصه لهمتمصتاهر نا-1 )عم خصمتعتمومن 
کصمحصتلعی عسفهه ما فصمتتومی حصیماتاتنومصمه تملصید 0۲)60حفصه مه واصمصصتلعو ملنطه ممتتصمی صیت‌طتاتناوم 
6 وه ۱۵۵۷۵ ممتانلجمن مصیت‌هاتآنناوه مط) مق م۱ ۵ تماههتو ۵۵ رقم فممتاتلجمم مصیاتهاتاتناوه م0 رمل‌صنا ات0مرقصقع 
0 ۱۵) ماهاباصله مه 0۵7۵100۵۵0 قه7 061مص لهمتتمصصباط مه بطم‌تهعوع: فنطا ما عمممصها مممصتمو که موه 
اممصصتلعه 060هممواو ۵ صتامصصع م1 .فممتاتلجم صا‌واتاننمع صمه عصا طمقع۲ تعتر و صا ا0۲حقصق اصممتلمو 
0 0 100باصاتاونل مط) فص 0ع7«مطاه فماناوع م1 بعمله صنهیع عصمحصتمو طمقع ۵۲ 0ماهلناماوه مه صمتحاجهعصمم 


۵21 روصتهعصتوم مماه ۷۲ ۵۶ مممانهمونا رممااگ ‏ ماه ۷۷ ۵۶ ماو ۱۷۲86 مه ۳۲۵1۲۵990۲ ملقاهموو۸ -2 20 1 
۰ ,1۳21 رطفطفصمصمک رنتفهب تون 

رصم تمصع ۵ صمتل‌هطامطعو :اتف تمه م مت ممصومتزمن -۶) 

10]: 6 


بررسی عددی عوامل موّثر بر توزیع غلظت رسوب معلق نامتعادل در رودخانه‌های طبیعی ۲۵۳ 


اصعصصتلعو ۵۶ 0۷۵۵ عظ) عانط/ دصمناقاه: باتهم(ع۷ 1211 عطا م) 1)1۷6عصعو باصم‌تصونه ام و1 جمقع بع۲۷ معط هعملة 1020 
۸ ۶ مومع فط1. مصمتامتامهعصمع )۲ممقصهتا اصمحصتلعو 060همم‌فتای فط) 260160 همتامنامع )225001 
ص واصمصتلعع ۶ حملاحتاجهعجمع فطا صقط . معطهنط 11.720 موله مه مصمتتقصمی مصیت‌ماتانتانعصمظ. رم واصمحصتلعو 

۰ صیات مات آنناوم 


0 ,20101 اوه لمآ بحمتتوم صی‌طتاننومعصما رطقمممممک هت ممفهتقطان :6۲۱۷۵0۵۵8 
اصمصصنلع؟ 


